Comment on "Discrepancies in the resonance-fluorescence spectrum calculated with two methods" 
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There are two alternative methods used in literature to calculate the incoherent part of the spectrum of light 
scattered by an atomic system. In the first, one calculates the spectrum of the total light scattered by the system 
and obtains the incoherent part by subtracting the coherent part. In the second method, one introduces the 
fluctuation operators and obtains the incoherent part of the spectrum by taking the Fourier transform of the two- 
time correlation function of the fluctuation operators. These two methods have been recognized for years as two 
completely equivalent for evaluating the incoherent part of the spectrum. In a recent paper Qing Xu et al. [Phys. 
Rev. A 78, 013407 (2008)] have shown that there are discrepancies between the incoherent parts of the stationary 
spectrum of a three-level A-type system calculated with these two methods. The predicted discrepancies can be 
severe that over a wide rage of the Rabi frequencies and atomic decay rates, the spectrum calculated with the 
variance method can have negative values. This is obviously unphysical result since the fluorescence spectrum 
is a positively defined quantity. It represents the frequency distribution of light incoherently scattered by the 
atomic system. Therefore, the calculated spectrum should be positive for all frequencies independent of values 
of the Rabi frequencies and the damping rates. In this comment, we show that there are no discrepancies between 
these two methods. The equivalence of these two alternative methods leads to the same incoherent spectra that 
are positive for all frequencies independent of values of the parameters involved. The analytical analysis is 
supported by simple numerical calculations. 
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The incoherent part of the fluorescence spectrum of an 
atomic system driven by lasers fields is commonly studied 
with two, often called "limit" and "variance" methods. In the 
limit method, one calculates the total spectrum from which 
then the coherent part is subtracted to obtain the incoherent 
part alone. In the variance method, one define the fluctua- 
tion operators and directly obtains the incoherent part of the 
spectrum by a Fourier transformation of the two-time cor- 
relation function of the fluctuation operators. These meth- 
ods are generally recognized as completely equivalent and 
have been used for decades to calculate spectra of the sta- 
tionary electromagnetic field radiated by variety of atomic 
systems tfH 0, HT H IH Hb . The variance method 

has been particularly popular in the calculations of the in- 
coherent spectrum of non-stationary fields due to some dif- 
ficulties in the treatment of the coherent part of the spec- 
trum 

In a recent paper, Xu et al. [17] have calculated the inco- 
herent spectrum of the fluorescence field emitted by a three- 
level A system driven by two laser fields and claim that these 
two methods produce different results, and the method of the 
calculation of the spectrum involving the fluctuation opera- 
tors may even lead to negative values of the spectrum. This 
statement, however, is incorrect. In this comment we clar- 
ify this point, directly showing that both methods lead to the 
same incoherent part of the spectrum and the negative values 
of the spectrum simply arise from an error in the Xu et al. lfl7ll 
numerical computation. Specifically, the explanation of the 
source of the discrepancies given by Xu et al. is based on the 
general solution of the set of differential equations describing 
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the time evolution of the system. Therefore, it is necessary to 
show that there are no discrepancies between these alternative 
methods since similar discrepancies could be found in other 
systems described by a similar set of differential equations. 

As it is well known, the spectrum of the fluorescence light 
emitted by an atomic system is obtained by a Fourier trans- 
formation of the two-time correlation function of the emitted 
electromagnetic field that is usually expressed in terms of the 
two-time correlation function of the atomic operators. In the 
case of the A-type atom considered by Xu et al. lfl7ll . the spec- 
trum of the fluorescence light emitted on the atomic |3) — |1) 
transition can be written as 
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where 71 is the spontaneous emission rate of the transition, 
u(r) is a geometrical factor corresponding to the radiation 
pattern of the dipole moment and the expectation values of 
the atomic dipole operators are evaluated in the rotating frame 
oscillating with the frequency lo\ of the atomic transition in- 
volved. 

Usually the scattered field is composed of a coherent com- 
ponent, corresponding to a field elastically scattered by the 
source atoms, and an incoherent (noise) component, corre- 
sponding to a field produced by fluctuations of the atomic 
dipoles. Therefore, it is often distinguished between these two 
contributions to the scattered field by expressing an atomic 
dipole operator as the sum of its expectation value and its fluc- 
tuations so that 

°-<i(*) = + Ao-yO), (2) 

where (Aer^ (t)) =0 by definition. In terms of the fluctuation 
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operators the spectrum is of the form 
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We illustrate the equivalence of the two methods using three 
different arguments that correspond to different approaches to 
evaluate the incoherent spectra. 

The first and the simplest argument to prove the equivalence 
of the two methods is to show that the incoherent part of the 
spectrum calculated with the "variance" method can be trans- 
formed to a form identical to the incoherent part of the spec- 
trum calculated with the "limit" method. To show it, we sim- 
ply expand the variance spectrum into two terms 
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where the first term on the right-hand side corresponds to the 
total spectrum. The second term corresponds to light that is 
elastically scattered by the atom. It is easy to see. Since in 
the steady-state (031 (f + r)) and (013(f)) are independent of 
time, we can write 

S in (oj) = S T (v) - 7iu(f) |(o3i) s | 2 S(u - ui), (5) 

where (031)3 is the steady-state value of the average atomic 
dipole moment of the |3) — |1) transition. 

We observe from Eq. dD that the resulting difference is man- 
ifestly positive since the coherent part of the spectrum cannot 
overweight the total spectrum. Thus, the incoherent part of 
the spectrum calculated with the variance method is neces- 
sarily positive. It is always positive independent of a mathe- 
matical approach adopted to evaluate the two-time correlation 
function of the fluctuation operators. In addition, note that the 
coherent part contributes only to the central component of the 
spectrum. Therefore, any discrepancies between the spectra 
observed at the Rabi sidebands do not come from the proce- 
dure of subtracting the coherent part of the spectrum. 

As the second argument of the equivalence between these 
two methods, we consider equations of motion for the expec- 
tation values of the atomic operators and show that the solu- 
tions for the atomic two-time correlation functions appearing 
in Eqs. (QJ and ® lead to the same incoherent parts of the 
spectrum. We consider the general case of systems whose 
dynamics are described by a set of inhomogeneous linear dif- 
ferential equations 

^X(f) = QX(t) + R, (6) 
at 

where Q is a finite-dimensional n X n matrix of time inde- 
pendent coefficients, Rabi frequencies, damping rates and de- 
tunings, involved in the dynamics of a given system, X(t) is a 



column vector of the one-time correlation functions and R is 
a column vector of inhomogenous terms. The dimension n of 
the vector X is equal to the dimension of the Hilbert space of 
a given system. The general case includes, as a special case, 
the dynamics of A-type atom driven by two laser fields, con- 
sidered by Xu et al. IU7I1 for which n = 8. 

We first find the steady-state solution of Eq. (0, which is 
obtained by taking the limit of f — ► 00, or more directly by 
setting the left-hand side of Eq. (|6]l equal to zero. Thus 

X(oo) = -Q^R. (7) 

Next, we consider two-time correlation functions involved in 
the calculation of the spectrum of a given system. These are 
simply found by applying the quantum regression theorem to 
Eq. © which states that for t > the two-time correlation 
functions satisfy the same equations of motion as the one-time 
averages 01811 . Thus, we find that 

-£-Y(t,T)=Q,Y(t,T) + (X i (t))R, (8) 

where Y(t, r) is a column vector of components 
(Xi^X^t + r)), (Xi(t)X 2 (t + r)), . . .. 

Using the Laplace transform, we can transform Eq. © into 
a set of algebraic equations, which we can solve by the matrix 
inversion. Thus, applying the Laplace transform, we obtain 

Y(s) = (.si - Q)- 1 y(0) + [(si) {sI-QT^XiityR, (9) 

where Y(s) is a column vector composed of the Laplace trans- 
forms of the correlation functions, s is a complex (Laplace 
transform) parameter and I is the unit diagonal matrix. 
We can rewrite Eq. (O in the form 

Y(s) = (sI-Q)- 1 (y(0) + Q- 1 (X i (f))^) 

-( s I)- 1 Q- 1 (X,(f))i?, (10) 
which in the steady-state limit of f — > 00, takes the form 
Y{s) = ( S I-Q)- 1 Ay(0) + ( S I)- 1 (X l (oo))X(oo),(ll) 

where AY(0) = Y(0) - (Xi(oo))X(oo). 

As we see from Eq. (fTTT > that the Laplace transform of the 
two-time correlation functions is composed of two terms, one 
corresponding to contributions of terms with s / that re- 
flect the presence of incoherent scattering resonances, and the 
other contribution from the pole at s = that reflects the 
coherent scattering peak. This term is proportional to delta 
function centered at the driving field frequency since the in- 
verse Laplace transform of the term proportional to (si) -1 
gives the delta function type contribution. 

We now show that the first term on the right-hand side of 
Eq. (fTTT i is equal to that obtained for the correlation functions 
of the fluctuation operators, i.e. leads to the same correlation 
function that is involved in the variance method of the calcu- 
lation of the incoherent part of the spectrum. This will con- 
stitute that both methods of calculating the incoherent part of 
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the spectrum are completely equivalent. We define a column 
vector of average values of the fluctuation operators AX(t) 
whose the time evolution is determined by the homogeneous 
equation of motion 

j t (AX{t)) = Q(AX(t)), (12) 

from which we find, by applying the quantum regression the- 
orem, that the vector AY(t, t) of the two-time correlation 
functions satisfies an equation of motion 

^AY(t,r) = QAY(t,T), (13) 

where the components of the vector AY(t, r) are 
(AX i (t)AX 1 (t + r)>, (AXi(t)AX 2 (t + r)>, .... 
Taking the Laplace transform of Eq. H3i , we find 

AY(s) = (si- Qy 1 AY (0), (14) 

which is equal to the first term on the right-hand side of 
Eq. (fTTT i. Thus, the limit and variance methods both lead to 
the same solutions for the correlation functions involved in the 
definition of the incoherent part of the spectrum. This state- 
ment is true for an arbitrary system whose the time evolution 
is described by the set of differential equations (|6). 

Third and the last argument illustrating that there are no dis- 
crepancies between the two methods is the numerical analysis. 
Figure Q] shows the incoherent part of the fluorescence spec- 
trum S3i(uj) calculated numerically with the variance method 
for the same parameters as in Fig. 2(b) of Ref. ifnll . Note that 
this is the figure where the most prominent negative values of 
the spectrum have been predicted. The negative values of the 
spectrum result from dispersive rather than Lorentzian-type 
structures at the Rabi sidebands. 




FIG. 1: The incoherent part of the fluorescence spectrum Ssi(uj) 
calculated with the variance method for the same parameters as 
in Fig. 2(b) of Ref. fl7tl . The spectrum is positive for all frequen- 
cies and exhibits Lorentzian-type structures at the Rabi sidebands. 



We plot the spectrum by taking the Laplace transform of 
Eq. (fT3l and solving numerically the resulting set of algebraic 
equations by matrix inversion 111911 . Evidently, the spectrum is 
positive for all frequencies and exhibits Lorentzian-type struc- 
tures at the Rabi sidebands, that is the incoherent spectrum is 
the same as calculated with the limit method, shown as a dot- 
ted line in Fig 2(b) of Ref. 0. 

In summary, we have presented three simple arguments 
against the statement of Xu et al. Il7ll that the incoherent part 
of the fluorescence spectrum of a three-level A-system calcu- 
lated with the variance method can have negative values. The 
arguments clearly show that the incoherent part of the spec- 
trum calculated with the variance method is positive for all 
frequencies independent of the system considered and values 
of the parameters involved. 
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